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clidean and hyperbolic spaces and are simultaneously approached starting 
from the Lie algebras 50^ (A + 1), which include a parametric dependence on the cur- 
vature. General expressions for the Hamiltonian and its integrals of motion are given 
in terms of intrinsic geodesic coordinate systems. Each Lie algebra generator gives rise 
to an integral of motion, so that a set of A(A + l)/2 integrals is obtained. Further- 
more, 2A — 1 functionally independent ones are identified which, in turn, shows that 
the well known maximal superintegr ability of the Smorodinsky-Winternitz system on 
is preserved when curvature arises. On both and H^, the resulting system can 
be interpreted as a superposition of an "actual" oscillator and A "ideal" oscillators 
(for the sphere, these are alike the actual ones), which can also be understood as A 
"centrifugal terms" ; this is the form seen in the Euclidean limiting case. 
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1. Introduction 

Superintegrable systems on the two- and three-dimensional (3D) Euchdean spaces 
have been classified in [6, 25], and also extended to the 2D and 3D spheres [11] as well 
as to the 2D hyperbolic plane [19, 20] and 3D hyperbolic space [12]. Recent classifica- 
tions of superintegrable systems for these 2D Riemannian spaces have been presented 
in [18, 21, 24, 27]. In particular, in the 2D sphere there are two (maximally) 
superintegrable families: the harmonic oscillator and the Kepler-Coulomb potential, 
both of them with some "additional" terms. Let us consider §^ as embedded through 
si + si + S2 = 1 in an ambient space = {sq, si, S2); we set the geodesic polar 
coordinates {r,6) such that sq — cosr, si = sin r cos ^, S2 — sin r sin ^. Following the 
notation and results given in [27], the first classical superintegrable family on §^ is 
given by 



(1.1) 



^ho — Po ( ^ _2 ^ ) + ^ + ^ 



Sq J Si S2 



= 13q tan r -|- — -|- 



sin^ r cos^ d sin^ r sin^ 9 



where (5i are real constants. The first term, tan^ r, is the spherical oscillator or Higgs 
potential [16, 26]. Under contraction to the Euclidean plane, this reduces to the usual 
harmonic oscillator, r^, while the two remaining terms give rise to two "centrifugal 
barriers" . However, very recently, all the three terms in Who have been interpreted as 
spherical oscillators with different centers [29]. 

On the other hand, the second superintegrable family on §^ turns out to be 

(^12) V«i + «2 sl^Jsl + sl si 

-R ^ ^R I 
— A^O 7 r Pi T-oT; + 



tan r sin r sin 6 sin r sin 6 

where 1/tanr is the "spherical" Kepler-Coulomb potential, first studied by Schro- 
dinger [30]; the two potentials tan^r and 1/ tanr are mutually related [22, 23]. 

The aim of this contribution is to study the maximal superintegrability of the 
generalization of the potential (1.1) on the A^D spaces S"^, and M.^ from a group 
theoretical standpoint. This family, depending on the curvature as a parameter, 
includes for k = the well known maximally superintegrable Euclidean Smorodinsky- 
Winternitz (SW) system [7-10], the Hamiltonian of which reads 

(1-3) n^\Y,(pl + 2p,^- + ''4)^ 

^ i=i V / 
where g| = is the harmonic oscillator potential and each l/g| is a "centrifugal 
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term". Two sets of integrals of motion for Ti. are given by (i < j; i,j — l,...,N): 
Ioi = P^ + 2Poq'^ + 2Pi\, with Pi=Pi, 

^ij = Jfj + + 2/3j ^, with Jij = QiPj - QjPi. 

The first set comes from the separabihty of the Hamiltonian 2H = ^j/oi, while 
the second one has a term quadratic in momenta (the square of the components of 
the angular momentum tensor) plus some additional terms depending on coordinates. 
The functions Pi, Jij close in phase space the commutation relations of the Euclidean 
algebra iso(A'") with respect to the canonical Lie-Poisson bracket: 

The structure of this contribution is as follows. The next section contains the 
details required (related to their maximal groups of isometries) on the three A^D 
classical Riemannian spaces with constant curvature k. Geodesic motion is obtained 
in section 3 starting from the metric on these spaces; this is the kinetic energy term to 
which possible potentials can be added. The generalization of the SW family to these 
spaces (any k) is performed in section 4. General expressions for the Hamiltonian 
and its integrals of motion are explicitly given both in terms of the (Weierstrass) 
coordinates in an ambient linear space M^+^ as well as by means of two sets of 
intrinsic coordinates, the non-zero curvature versions of the Euclidean Cartesian and 
polar ones. Moreover, — 1 integrals, including the Hamiltonian, are shown to be 
functionally independent, thus proving that the SW system on curved spaces is also 
maximally superintcgrablc [4]. As an example, we apply in the section 5 the general 
expressions to the = 4 case. 



2. The sphere Euchdean E^, and hyperbolic spaces 

Let S0k{N + 1) be the Lie algebra of the motion group SOk.{N + 1) on a generic 
A'"D real Riemannian space with constant curvature k, denoted 5^]. In the basis 
{ Joi = Pi, Jij} = 1, . . . ,N; i < j), the non-vanishing commutation relations of 
50k{N + 1) are given by 

^2 [Jijj Jik] — Jjki \'^iji '^jk] — Jikj [Jikj Jjk] — Jiji 

[Jij J -fi] — Pjj ['^iji Pj] — [Pii Pj] — I^Jiji 

with i < j < k. The Lie algebra 50k (iV + 1) is isomorphic to either so (A^ + 1) for k > 0, 
iso(A'") for k; = 0, or so{N, 1) for k < 0. Notice that by scaling the generators Pi, any 
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value of K can always be reduced to either +1, or —1. The involutive automorphism 
defined by 

(2.2) B : Jij Jij, Pi -Pi, i,j = l,...,N, 
provides the following Cartan decomposition of so^iN +1): 

(2.3) 50«(7V + l) = f)®p, i) = {J,^)=5o{N), p = {Pi). 

The generators invariant under span the subalgebra f) of the Lie subgroup H ~ 
SO{N), so that S^-^ = SO^{N + 1)/S0{N) is a family, parametrized by k, of iVD 
symmetric homogeneous rank-one spaces [13]. The generators Jij leave an origin 
point O invariant, thus acting as rotations around O, while the remaining Pi generate 
translations that move O along basic geodesies li orthogonal at O. The space S'^j 
comprises the three classical Riemannian spaces: 

K > 0, Sphere, 3(1^ = SO{N + 1)/S0{N) = 

K = 0, Euclidean, 5g = IS0iN)/S0iN) = E^; 
K<0, Hyperbolic, = SO{N, l)/SO{N) = H^. 

The curvature n can also be interpreted as a graded contraction parameter coming 
from the Z2-grading of 50k.{N + 1) determined by G; the value k = corresponds to 
the contraction E>^ — > E'^ <— H'^ around a point (the origin O) of the spaces. 

The Killing-Cartan form of so«(Ar + 1), g^'^{Jij, Jki) = Trace(ad J^j ■ ad J^z), is 
diagonal with the following non-zero elements 

(2.4) g'^^iPi, Pi) = -2(iV - 1)k, g^'^iJij, Jij) = -2{N - 1); 
hence the restriction of g^'-^ to the subspace p, can be written as: 

(2.5) g''^\= -2{N - l)Kg\^, g\^{Pi, Pj) = Si^, 

where 5f|^ is to be considered as the metric in the tangent space p at C, represented 
by the N x N matrix diag (-I-, +,..., +). This metric can be translated to all points 
of /Sj^j by the group action. Even if g^^ vanishes in the Euclidean case, the choice 

g oc g^'^ / K. ensures a non-degenerate metric in jS^j Vk. 
2.1. Vector model and Weierstrass coordinates 

The vector representation of 50^{N -\- 1) is given by (A^ + 1) x (A^ + 1) real matrices 
fulfilling (2.1): 

(2.6) Pi = KCQi -\- 6jO) Jij — ~l~ Cjij 
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where is the matrix with a single non-zero entry 1, at row i and column j. Their 
exponential gives rise to the following one-parametric subgroups of SOk.{N +1): 



AT 



exp{a;Pi} = ^ e^s + C«(a;)eoo + C«(a;)eii - KS«(a;)eoi + Sf^{x)e. 
ex-p{xJij } = Bss + cos X en + cos x ejj — sin x Cij + sin x e 



iO 



s=0;s7^i,ji' 



= I -h sin X + Jfj{l — cos x) , 

where I is the {N + 1) x {N + 1) identity matrix. The curvature-dependent cosine 
Ci^{x) and sine S^{x) functions are defined by [3]: 



x 



21 



1=0 



(2.8) 




X 



21+1 



Sm^/KX, 



/=0 



{2l + l)\ 



X, 



K > 0; 
At = 0; 



sinh V— K^c, K < 0. 



From them, we define the "versed sine" or versine V i^{x) and the tangent T^(x): 
(2.9) 



V.(x) = -(1- C.(x)), T^{x)^^^ 



These K-trigonometric functions coincide with the usual elliptic (circular) and hyper- 
bolic ones for k = +1 and k = —1, respectively. The "flat" case with k = gives the 
parabolic functions: Cq{x) = 1, Sq{x) = x and Yq{x) = x^/2. In this sense, they can 
be interpreted as the "curvature n-deformed versions" of 1, x, and x'^. All the known 
trigonometric relations (necessary in the further computations) also extend for these 
K-functions [14, 15] such as, for instance, 



— Ci^{x) = -K S^{x), — 
dx dx 



Si^{x) = Ci^{x). 



In the vector representation (2.6), any generator X satisfies the relation 



JV 



(2.10) 



X'^I^ + I^X = 0, Ik = eoo + i^'^eii= diag(l,K, ...,k), 



i=l 
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(X^ is the transpose matrix of X), so that any element R e SO^iN + 1) verifies 
R^I^R = Ik- In this way, SOk.{N + 1) can be seen as a group of hnear isometries 
of the bihnear form 1^ acting on R^"'"^ = (sq, si, . . . , sjv) via matrix multiphcation. 
The action of SOk{N + 1) on is hnear but not transitive, since it conserves the 

quadratic form 

N 
i=l 

provided by and H ~ SO{N) = (Jij) (2.3) is the isotopy subgroup of the point 
O = (1,0, ...,0) G R-^"'"^ which will be taken as the origin in the space The 

action becomes transitive if we restrict to the orbit in M^+^ of the point O, which is 
contained in the "sphere" E: 



N 



(2.11) 



S = sl + K^^. 



1=1 



which reproduces the whole hypersphere, two hyperplanes and a two-sheeted hyper- 
boloid for k >, =, < 0, respectively. This orbit is identified with the space S'^p and 
{so, si, . . . , sn), fulfilling the "sphere" constraint (2.11), are called Weierstrass coor- 
dinates. In terms of these, the metric on 5'^^^ comes from the flat ambient metric in 
■[^iv-i-i (jiyi(jg(j ^ restricted to S: 



N 



(2.12) 



dc7^ = - I dso + K ^ dSi 



i=l 



N 



1=1 



« Eili dsl 



which reduces to the Euclidean one for k = 0. A differential realization of the gener- 
ators as first-order vector fields in E-'^+-'^ with di = d/dsi is (see (2.6)): 



(2.13) 



K Sido - sodi 



Sjdi Sidj. 



2.2. Geodesic coordinate systems 

Let us consider in the vector model the origin O = (1, 0, . . . , 0) and a generic point 
Q e S'j^j with Weierstrass coordinates s = {sq, si, . . . , sn) £ M''^"'"^. Starting from O, 
the point Q can be reached in different ways through the action of N one-parametric 
subgroups (2.7), that is, by means of motions in the space S^y In this way, we intro- 
duce coordinates which are intrinsic quantities to the space S^^^ itself, the associated 
Weierstrass coordinates automatically fulfilling the condition (2.11). We shall consider 
two possibilities. 
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2.2.1. Geodesic parallel coordinates a = (ai, . . . , oat) . 

We move the origin O by using the N translations subgroups as 

s(a) = exp(aiPi) exp(a2P2) • • •exp(aAr-PAr) O 

( C^(ai) C«(a2) 0^(0,3) ■ ■ ■ C^(aiv) \ 



(2.14) 



/ so \ 












SN-1 









S«;(ai) C^(a2) 0^(03) . . . C^{aN] 
8^(02) C«(a3) . . . C«(aiv) 

Sfj(aAr_i) Cfj(aiv) 

V SK(aAr) 



/ 



Each coordinate is associated to the generator Pi and has dimensions of length. If 
/i, /2, . . . , /at are the (oriented) basic geodesies in S'j^j orthogonal at (9, then the first 
coordinate of a point Q is the geodesic distance ai between O and Qi (the orthogonal 
projection of Q on Zi), measured along Zi. The second coordinate is the distance 
02 between Qi and another point Q2, measured along a geodesic I2 orthogonal to li 
through Qi and parallel to I2 (in the sense of parallel transport) and so on, up to 
reaching Q [15]. This is depicted in figure 1 for and tf. By introducing (2.14) in 
(2.12), we obtain the metric: 



(2.15) 



N-l / N 

d^'= E n Glial) ]da^ + da 



From it we may compute the Levi-Civita connection, whose non-zero Christoffel sym- 
bols are (i = 1, . . . , iV - 1; i + 1 < j < N): 



(2.16) 



l=i+l 



2.2.2. Geodesic polar coordinates 6 = {r,92, ■ ■ ■ ,0n). 

In this case, we move O by using the N—1 rotations Jj j+i as well as the first translation 
Joi = Pi: 

s{0) = exp{eNJN-iN) ■ . .exp(6'2Ji2) exp(rPi) O 
/ So \ / C^(r) 
si S«(r)cos^2 
(2.17) S2 S«(r) sin 6*2 cos 6*3 



V Sjv / 



Si.-{r) sin 6*2 •• • sin 6* at- 1 cos 6* at 
\ S„(r) sin ^2 • • .sin^Ar-i sin^Ar / 
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Figure 1. Geodesic parallel (01,02) and polar (r, ^2) coordinates in 
and IP. 

The "radial" coordinate r associated to Pi has dimensions of length and is the distance 
between O and Q measured along the geodesic I joining both points. The remaining 
$2, ■ ■ - jOn are ordinary angles parametrizing the orientation of I with respect to the 
reference flag at O spanned by {/i, /2, • • • , ^i-i} (see figure 1 for = 2). On the sphere 
S'^ of radius R and curvature k, — all the usual spherical coordinates are angles, 

and these differ from the geodesic polar ones only in the first coordinate [28], which 
is commonly taken as the dimensionless quantity di = r/R. While these conventional 
spherical coordinates require an explicit contraction to the Euclidean ones, as done 
for instance in [17], our choice works for any value of k and when k = we recover 
directly, without any limiting procedure, the polar (and Cartesian) coordinates in E-^. 
In polar coordinates the metric (2.12) turns out to be 

(2.18) da^ = dr^ + (r) (dej + ^ ( H ^i) 

The components of the Riemann and Ricci tensors can be computed in either 
geodesic coordinate system; the scalar and the sectional curvatures of S^^^^ are R — 
N{N — 1)k and K — k,, both constant, respectively. When k — 0, the expressions 
(2.15) and (2.18) give the Euclidean metric 

(2.19) da^ = ^ da\ = dr^ + d6'^ + ^ ( sin^ j dOl 
and all the Christoffel symbols (2.16) vanish in parallel coordinates. 
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3. Geodesic motion on S^^ and phase space realization of so^(iV+l) 

If we now adopt a dynamical viewpoint, the expressions of the metric in S'j^] (2.15) and 
(2.18) provide the kinetic energy T of a particle in terms of the velocities, generically 
denoted as q (either d or 9) in the coordinate systems (2.14) and (2.17): 

7V-1 . N X / N .1-1 x\ 

(3.1) 2r = ^ n ^"Mj al + d% = r'+ S^(r) + E 11 f H ' 

which is the Lagrangian £ = T of the geodesic motion on S^y The canonical mo- 
menta, p = dC/dq, (denoted p, n for q — a, 9) are given by: 



Pi={ W'^'i{ai)\ai, i = l,...,N -1; p^ = dN, 
\i=i+i J 

7ri = f, 7r2=S^(r)^2, tt,- = S^(r) (j] sin^ j = 3, . . . , iV. 



d=2 



By introducing these momenta in (3.1), we obtain the free Hamiltonian 7f = T in the 
phase space of motions in S'j^j expressed in either "parallel" (a,p) or "polar" (6',7r) 
canonical coordinates and momenta: 



N-l 2 1 / ^ 2 \ 

(3-3) 2r=E ^ _ +^i^--^^+s4^^i + E ^^-^^ 2j - 

On the other hand, the Lie generators are expressed in terms of Weierstrass coor- 
dinates as: 

(3.4) Pi{s{q),s{q,p)) = SQSi- SiSo, Jij{s{q), s{q,p)) = SiSj - SjSi, 

and thus we get an A'"-particle realization of so^{N + 1) in the phase space simply 
by rewriting everything either in terms of {a,p) or (^,7r). The time derivatives of s 
are obtained from either (2.14) in terms of parallel coordinates and velocities (a, d) 
(i = l,...,iV-l): 

N N 

Tn=l 1=1 
N / N \ 

Yl ^Kictm) I CLi - KT^(ai) E T^Kiai)ai J , sn = C«(ajv)dAr, 

Tn=i \ l=i+l / 



So = 



(3.5) 
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or from (2.17) in terms of polar coordinates and velocities (6', 6) {j = 2, . . . ,N — 1): 



so = -KS^{r)r, si = S«(r)sin6'2 7^7— — — - 0^ 

i+i ( r ^ e 

(3.6) h = S.(r) n sin^^ T.(r) tan^,,., + g tan^.tln^,^, " '^'^^ ) ' 

i.= S.(r)n-^.. + 

We now introduce the parametrizations (2.14), (2.17), velocities (3.5), (3.6) as well as 
the momenta (3.2) in (3.4), obtaining a phase space realization of the generators of 
S0k{N + 1) given in "parallel canonical" coordinates hy = 1, N): 

Pi=il CM CMPi + «S,(a,) S.(a,)n^p^^ 



m I 



(3.7) 4- = S.(a.) CJa,) n ^Mv, " ^^^^^^^^,Vi 

I c (' ^ c (' ^ c ( \ nm=i+i CK(o.m) 

+ «S«(ai)S«(aj) 2^ S«(as)— — — Ps, 

and in geodesic polar coordinates and momenta (6*, tt) by {i,j = 1, . . . , N — 1): 

~ nSsin^fc UmLs^^^^rn cos dsTTs TT^+i 



tan^i+i ^ T«(r)tan^,+inL2sin^/ T«(r)nSsm^r 

T . /I /I TT • /I cos^i+i sin^^-i_i 
(3.8) = sme/i+i cose/j+i II sme/feTTi+i —■ ^ tt^+i 

+ COS 6^i+l COS 6/^+1 2^ -— TTs, 

T n n n , /) 1 1^^^ sm 6*^ COS 
JiN = sindi+i sm 6'fc TTi+i + COS 6'i+i 2^ — = — tt^. 

k=i+i s=i+i i-h=i+i^^^^i 

Then the following results follow [4] . 

Proposition 3.1. Both sets of generators (3.7) and (3.8) fulfil the commutation rules 
(2.1) of 50k.{N + 1) with respect to the canonical Poisson bracket (1.5). 
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Proposition 3.2. Any generator (3.7) and (3.8) Poisson- commutes with the kinetic 
energy T (3.3). 

The second statement is straightforward as T is related with a phase space real- 
ization of the second-order Casimir C of sOk{N + 1) through 

AT N 

(3.9) 2r = c^Y.P'^^Yl 4- 

In fact, the geodesic motion is maximally superintegrable and its integrals of motion 
come from any function of a phase space realization of the Lie generators. 



4. Smorodinsky— Winternitz system on curved spaces 

The problem now is to find which potentials U{q) can be added to T in such a manner 
that the new Hamiltonian Ti = T ^lA preserves the maximal superintegrability. This 
requires to add "some" terms to "some" functions of the generators in order to ensure 
their involutivity with respect to Ti. 



4.1. Potential 

As far as the 2D potential (1.1) is concerned, the generalization to the space S^^^ is 
straightforward and reads 

(4.1) W(.(g)) = /?o^^ + E3, 

which in geodesic parallel and polar coordinates turns out to be 

^^"^^ Slir) [cosH, + ^ cos^^.+inUsin^^. ULsln^O^^ 

This potential, which coincides with (1.1) for k = 1, has been interpreted on as 
three Higgs spherical oscillators, whose "centers" are placed at the three vertices of a 
sphere's octant [29] (recall that the Higgs oscillator has two antipodal "centers"); this 
interpretation can straightforwardly be extended to the sphere with k, = 1/i?^ > 0. 
Let Oi be the points placed along the basic geodesies k and a quadrant apart from O 
(each pair taken from 0,0i are mutually separated a quadrant distance n/{2^/K) = 
i?7r/2 on S^). If we set k = 1 and denote r, the distances between any generic point 
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Figure 2. Distances involved in the SW potential on and H^. 

Q and (9, Oi along the joining geodesies, we find sq = cos r, Si = cos (i = 1, . . . , A?"), 
and (4.2) can be rewritten as 

N n N N 

(4.3) W = /?otan2r + J]— 9— =/3otanV + ^/?itanVi + J^A, 

COS 

i=i * i=i i=i 

which can clearly be recognized under this form as the joint potential due to a set of 
A'^ + 1 spherical oscillators whose centers are at the N + 1 points O, Oi. Alternatively, 
if Xi = it/2 — Ti, then each term in (4.3) can be described, as usual, as the 
spherical "centrifugal" barriers Pi/ sin^Xj. Under the contraction k = 0, — > E^, 
the Higgs-term gives rise to the "flat" harmonic oscillator /5or^ = Po Yli while the 
N remaining oscillators (whose centers would be now "at infinity" ) leave the Euclidean 
"centrifugal" barriers Pi/xf = Pi/aj as their imprints. 

On the hyperbolic space the potential (4.2) can similarly be interpreted. Let 
Xi be the distance between the generic point Q and the totally geodesic codimension 
one submanifold through O orthogonal to the geodesic Zj. When k = — 1 we find that 
So = coshr. Si = sinha^i (i = 1, . . . , N), so the potential can be written as: 



(4.4) U = Po tanh^ r + V ^1 = Po tanh^ r + V — % V A. 

^ sinh Xi ^ tanh Xi ^ 

1=1 * 1=1 * 1=1 

The first term is an "actual" hyperbolic oscillator with center at O (whose potential 
is bounded in the whole hyperbolic space), while each of the N remaining terms 
can be interpreted, as done conventionally, as some kind of hyperbolic "centrifugal" 
potentials, but also as "ideal" hyperbolic oscillators, whose centers would be beyond 
infinity, that is, in the exterior region of the hyperbolic space. 



13 



4.2. Integrals of motion 

By taking into account the results given in [27] for the integrals of motion of (1.1), let 
us consider the following functions, with quadratic dependence on momenta, given in 
the ambient space R-^+-'^ by (z < j, i,j = 0,l,...,N): 

(4.5) lij = {siSj -SjSif + 2(3i -| + 2(3j . 

Therefore we have A^(A^ + 1)/2 phase space functions coming from the Lie generators, 
which in the geodesic parallel canonical coordinates (3.7) turn out to be 

T - p2 I g.^ , UU Glial) 

(4.6) '-^ 



The same quantities read in the geodesic polar canonical coordinates (3.8): 
2 , T2(r)nil2sin^^/ , tan2^,+i 



(4.7) 



^ ^ tan^^... ^ T^J.) sin^ 0, ' 

^ 1 

r=2 T2(r) 11^=2 sm'^i 

r2 , cos2 ^,+1 nL,+i sin^ 2/3, COS^ ^,+1 



" ^ cos^ + eos^ 9,^, nL.+i sin^ 9^ ' 

~o n^,-^i sin^ 9i cos^ ^,+1 



1- ZPJV — -jT^^ 

cos2 9,+^ Yll.^^ sin^ 9i 



A first property for these quantities is given by: 

Proposition 4.1. The N{N + l)/2 functions given by either (4-6) or (4-V fulfil 
{lij, Iim} — whenever all four indices i < j;l < m, are different. 

From now on we consider the Hamiltonian H = T + U with T and U given in (3.3) 
and (4.2). Then we find that: 

Proposition 4.2. The N{N + l)/2 functions either (4-6) or {4-'^} ore integrals of the 
motion for the Hamiltonian, that is, {Iij,7i} — Wij. 
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We remark that the property analogous to (3.9) is now given by 

N N N 

(4.8) 2n = J2Io^ + ^^ Yl ^^j+2«I]A- 

i=l i,i=l «=1 

Notice also that when k = 0, the Hamiltonian expressed in parallel coordinates and the 
integrals (4.6) with generators (3.7) directly reduce to the flat SW system characterized 
by (1.3) and (1.4). 

4.3. M£Lximal superintegrability 

The last step is firstly to identify, within the initial set of N{N + l)/2 integrals of the 
motion together with the Hamiltonian, N which are functionally independent and in 
involution in order to prove the complete integrability of 7Y, and secondly to find out 
how many are functionally independent thus analyzing its superintegrability. 

Let us choose two subsets of — 1 integrals Q^''^ and Q^^) {I — 2, ... , N) obtained 
starting from the integrals lij associated to the rotation generators as: 

I N 

(4.9) Q» = ^(0 - E ^^i' 

i,j = l i,j=N-l+l 

which share the element Q^^^ = Q(^n)- The generators that determine the quadratic 
terms in the momenta in the first set Q*^^-* are associated to a chain of orthogonal 
subalgebras within {) = so{N) = (Jij) (2.3) starting "upwards" from so(2) = (J12); 
likewise for but starting "backwards" from so{2) — (Jn-in)'- 

Q(') c...cQ^') c...cQ^^) Q^N)^---^ Q{i)^ •••3 Q(2) 

S0{2) C . . .C S0{1) C . . .C S0{N) S0{N) D . . . D S0{1) D . . . D 5o(2) 

We remark that the SW Hamiltonian on E-^, Ti,] can also be constructed 
through a coalgebra approach [5] by means of N copies of sp(2,]R) ~ s[(2,]R) [1]. In 
this fiat case, each integral of motion Q(^) | „ in Cartesian coordinates coincides with 
an A^-particle phase space realization of the k-th order left-coproduct of the Casimir 
of 5l(2,R); similarly for using the right-coproduct (these are called left- and 

right-integrals, respectively). Such coalgebra method shows that the N functions 
{Q*'''* |^_o' ^L-oJ^ are mutually in involution (the right-set as well), and the 2A'" — 2 
functions {Q^''^ L-o' L-o' ^L-oJ" functionally independent [2]. For the general 
case with arbitrary curvature n, we have the following result. 
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Theorem 4.3. (i) The N functions 

are mutually in involution. The same property holds for the second set 

{Q{N)i • • • 5 Q(3)) Q{2),'H}. 

(a) The 2N — 1 functions 

Q^^\ Q(^) = Q(^), . . . , Q(3), Q(2), loi, n}, 

where i is fixed, are functionally independent, thus 7i is a maximally superintegrable 
Hamiltonian. 

Notice that if k 7^ 0, each Q^''^ (or can be seen as a smooth function on the 
curvature: 

The same holds if one introduces a quantum deformation of the Euchdean SW system 
in terms of a deformation parameter z [1, 2]: superintegrability is preserved and the 
deformed integrals Qz^ (or Q^^)) can be written as power series of z: 

Then, Q^'-'>\^^^= Qz^\^^q are the integrals of the SW system on E^. This suggests 
some kind of duality between quantum deformations (z) and curvature (k) [3] . 



5. The N = 4 case 

In what follows we illustrate the general expressions obtained in the previous sections 
by applying them to the case N = 4 in. terms of geodesic parallel coordinates and 
momenta. The SW Hamiltonian on S^, E'^ and reads 



(5.1) 



a , rr.2r ^ T2(a2) T^lM T^liaA) 

+ T;(ai) + 



C^Jai) Ql{ar)Gi{a2) Cliai) Clia-i) Gliaz) 



S^.(«i)C^(a2)C2,(a3)C2(a4) (a^) C^Jas) (a^ 



S2(a3)C2(a4) S2(a4)' 
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The phase space reahzation (3.7) for 30^(5) is given by four translations 
Pi =pi, P2 = Gi^{ai)p2 + KS^{ai)T^{a2)pi, 
P3 = C«(ai) C«;(a2)p3 + «T«(a3) ( \ pi + C«(ai) S«(a2)p2 



C«(a2) 

(5.2) ^ / S^(ai) 

P4 = C«(ai) C«;(a2) 0^(03)^4 + KT^(a4) ^ , ^ , — rpi 

V 0^(02) 0^(03) 

, C.(ai) S.(a2) ^ ^^^^^^ ^^^^^^ S.(a3)p3^ , 



C,(a3) 

together with six rotation generators: 

J\2 = S«(ai)p2 - C«(ai)T«(a2)pi, 

^13 = S«(ai) C^(a2)p3 - ^""1"^^ T«,(a3)pi + k S«(ai) 8^(02) T^(a3)p2, 
-^14 = S«(ai) 0^(02) C«(a3)p4 - ^ . "i, — r TKia4)pi 
^^■^^ +kSM^M f|4^P2+ C«(a2)S«(a3)p3l , 



J23 = S«(a2)p3 - C«(a2) T^(a3)p2, 

^24 = 8^(02) C«(a3)p4 - ^ T«(a4)p2 + « 8^(02) 8^(03) T^(a4)p3, 

= 8« (03)^4 - C^{as)T^{aA)p3. 
Prom these generators we construct ten integrals of motion (4.6): 

Pi P2 P3 Pa Iqi I 02 Iqz Iqa 



J12 J13 Jl4 =^ I12 Il3 -^14 

J23 J2A I23 I2A 

J3A I3A 

Hence from (5.2) we obtain four "translation-like" integrals given by 

>2 I o/a rr2' X ~- 1 



/oi = Pi +2/3oT^(ai)+2/3i 



(5.4) 



/o2-P2+2/3o-^j^+2/32^^, 

T - 52 , ma3) , C^.(ai)C^(a2) 

Jos — -f3 + ^POT^TT TT^TT T + ^^3 I^^TT ^ : 

C^(ai)C2(a2) T2(a3) 

'04 = -f^4 ' 



p2 , T2(a4) C2(aOC2(a2)C2(a3) 

+ C2 (aO C2 (a2) C2 (a3) + (a4) 
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together with six "rotation-hke" ones coming from (5.3): 

Jia — -'la + — 77^27 — ^ + — ^ : 



(5.5) 



, _j2 , T^.M , S^Jai)C^,(a2)C^(a3) 

Yi4 - Jx4 + 2/^1 g,^^^^ ^^^^ + 2/^4 ^^^^ 



T^(«3) , S^,(a2) 
S^(a2) T^(a3) 



, _ 72 , T^.K) , S^(a2)C^(a3) 

^lMCl{ai) T2(a4) 

t2 nn '-'-'k(^4) ^kI'^s) 
/34 = J34 + 2/^3^1^ + 2/34^^. 

Thus, in this case, within the subset of "rotation-hke" integrals, we find three 
"upwards-integrals" Q'-'^\ Q'^^\ Q'^'^\ associated to so(2) C so(3) C so(4), that deter- 
mine a completely integrable Hamiltonian: 

g(2) = Q(3) ^ + /l3 Q(^) = Q(4) = /l2 + /l3 + hi 

+-^23 +-^23 + ^24 

+^34 

These three functions together with the two "backwards-integrals" Q(2),<5(3), associ- 
ated to so (2) C so (3), and one additional "translation-like" one, say /qi, characterize 
a maximally superintegrable Hamiltonian (5.1): 

Q{2) — h4 Q{3) — I23 + I24 loi 

+I34 

To end with, we would like to point out that a similar algebraic construction 
may likely be applied to the curved version of the "Kepler-Coulomb" superintegrable 
family (1.2). Furthermore, the consideration of a second contraction parameter, say 
K2, that enables to take into account indefinite metrics of Lorentzian signature [14, 
15], would allow one to obtain superintegrable systems on different spacetimes. The 
application of these maximally superintegrable systems in quantum mechanics also 
deserves a further study. 
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